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ADDENDUM TO USERS GUIDE TO VI ? ASA 

(Vibration and Instability of Plate Assemblies including Shear and Anisotropy) 

By 

Melvin S. Anderson, Katherine W. Hennessy 

and 

Walter L. Heard, Jr. 

Langley Research Center 


SUMMARY 


Extensions developed at Langley Research Center to the VIPASA computer 
program are described including a procedure for simple thermal stress analysis 
and options for graphical display of output. Input requirements for operation 
of the modified program are given in detail. 


INTRODUCTION 


The VIPASA program (References 1 and 2) for panel buckling and vibration 
analysis was originally made available to NASA by the United Kingdom. It was 
found to be a very useful, efficient and accurate analysis tool. Extensions 
developed at Langley for plotting, thermal stress, and orthotropic laminate 
property input have further enhanced its capability. The program is now 
generally available and the purpose of this paper is to document the changes 
and extensions so that this report and the original VIPASA users manual will 
serve as the users manual for the current version of the program. 


Extensions to the VIPASA Program 

A number of extensions to the VIPASA program which have been found to 
enhance its capability and efficiency have been incorporated at the Langley 
Research Center. These extensions are summarized below: 

1. Input of a uniform axial strain so that individual plate loads can 
be automatically calculated for all plate structures. 

2. Calculation of total axial load on the structure for each eigenvalue. 

3. Input of layer temperatures so that simple constant strain thermal 
buckling problems can be handled. 
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4. Input of individual orthotropic layer properties with internal 
calculation of plate stiffnesses. 

5. Calculation of geometry including consistency checks for sub 
structures and final structures. 

6. Plotting of undeformed shape and buckling or vibration mode shapes. 
Plotting of eigenvalue as a function of wavelength. 

7. An option to speed up the calculation process when only the minimum 
buckling load over a range of wavelengths is desired (It is only 
necessary to perform one iteration at a given wavelength if the eigen- 
value is higher than that determined for some previous wavelength.) 

8. Addition of a title card. 

The details of the thermal stress calculation and stiffness properties of 
layered wall configurations ar-i described briefly in Appendix A. Examples of 
plotting and geometry checks will be covered in the section entitled EXAMPLE 
PROBLEMS. 


Program Operation 

The modified program input is changed very little : om the original 
VIPASA input with the exception of layered property iripit if that option is 
used. Input for the original program will also execute the current program 
if a title card is added prior to the C cards for each problem number. 

The plotting and geometry calculation program is separate from VIPASA. 

The usual mode of operation is to run both programs in sequence. The complete 
input for the plotting program is generated in VIPASA and stored on an output 
file which becomes the input file for the second program. There are no other 
input requirements. 

Because plotting commands are generally tied to facilities, the different 
subroutines that execute plotting are described with comment cards in the 
source listing. There is enough information given so that available 
routines that perform the same function can be easily substituted. In addition, 
a more detailed description of each plotting subroutine is given in Appendix B 
so that definitions of all the variables may be obtained. The programs are 
written in FORTRAN and do not use any complicated instructions so that they 
should be capable of running on a variety of machines with little or no 
modification except for the plotting instructions. Experience at Langley 
with the FTN compiler on CDC hardware has been that OPT = 0 is required. 

The program will compile at OPT = 1 but will give incorrect results. 
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INPUT REQUIREMENTS 


The following section gives a description of the card input in the 
correct sequential order including the input for the extensions developed 
at Langley. The cards that remain unchanged will be dealt with rather 
briefly and the reader is referred to the original VIPASA users manual 
(Reference 1) for a more complete description when required. New cards or 
cards that have been changed are indicated by underlining the new input 
variables. 


Card Input 

The input is divided into several categories denoted by letters of the 
alphabet A through I with one or more cards in each category. 

GROUP A Sets up array dimensions and plotting options, 

NMK, NSR, NSC, IPAST, NSN, NSK, NOLAM, IPLOT, NSPAC, PAGE, AMP 
(717, 214, 2E10.4) 


NMK Must be even and no less than highest plate number or doubly 

connected substructure number and less than or equal to 120. 

Ml'sL be double this value with shear or if either D13 or D23?*0. 

NSR Must be even and no less than the largest number of nodes 

comprising a substructure. Must be no less than the highest 
numbered node of the final structure (twice this value for shear 
or anisotropy). Must be at least 4. 

NSC Must be greater than the largest diti. rence between connected 

node numbers and at least 3. 

IPAST Must be no less than the number of iterations required for one 

eigenvalue solution. 

NSN Must be no less than the total number of all the integers in 

group H and I cards. 

NSK Must be no less than the highest numbered singly connected 

substructure number minus 120 and less than 900. 

NOLAM Must be no less than the number of different values of X for 

which solutions are required. 

56*NMK + 16*NSR*NSC + 4*NSR + 2*IPAST + NSN + 18*NSK + 4*N0LAM + 9 < IDUM 

IDUM is currently set at 8009; see comment cards at beginning of VIPASA 


I PLOT 
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Controls geometry calculations and plotting 
0 No geometry calculations made and no plotting done 

VI Used for checking input. Geometry calculated and 

undeformed structure plotted. No eigenvalue calculations 
made. 

+2 Plots are made of undeformed structure and for each 

eigenvalue the buckling or vibration mode superimposed on 
the undeformed structure is plotted. 

+3 Plots are made of undeformed structure and for each 

eigenvalue, the buckling or vibration mode is plotted alone. 

A negative value of IPLOT causes additional printout giving all of the 

coordinates of the substructures and the final structure. 

If eigenvalues are calculated at more thajn one t value of wavelength A . 

a plot of FACTOR(Eigenvalue) versus A is also made for |IPL0T( > 1 

NSPAC Number of intervals between nodes where mode shape is calculated 

for plotting purposes. May be omitted and default value of 4 
will be used. Larqe values cause smoother mode shapes but no 
increase in accuracy. 

PAGE Controls size of plots. Overall dimensions of plots vary 

linearly with PAGE. If PAGE = 1., all plots will fit inside a 
6x9 rectangle (this may be changed in the program as indicated 
in comment cards by changing the definition of variables PX 
and PY). May be omitted and default value cf 1. will be used. 

AMP Maximum amplitude of buckling or vibration mode. May be omitted 

and default value of .5*PAGE will be used. 

GROUP B 

NOPROB (15) Number of problems. If more than one problem is desired, the 
A and B cards are not repeated for subsequent problems. 

GROUP C 

TITLE (I) (8A10) 

May contain any Hollerith text 

JBMAX, JBMIN, CONVER, FACLIM, FACLOW, SIGMA, ANV, MODE, LPLATE 

(215, F15.1 , 2F15.6, 2F10.4, 212) 

If JBMAX positive all eigenvalues with eigenvalue numbers 
between and including JBMIN and JBMAX will be found. If 
JBMAX negative all eigenvalues having a value between FACLOW 
and FACLIM will be found. (In this case JBMIN should be left 
blank) . 


JBMAX 

JBMIN 
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CONVER 

FACLIM 


FACL.0W 

SIGMA 

ANV 

MODE 

LPLATE 

EL, ALAMIN, 
EL 

ALAMIN, 

ALAMAX, 

RATLAM 


EC, ANUC, 
DENSEC 


Eigenvalues will be found to an accuracy of EIGENVALUE/CONVER 

Should be set to maxirr im expected eigenvalue if JBMAX is positive 
(does not have to be accurate but good guesses improve 
computational efficiency). Should be set to upper limit of 
desired range of eigenvalues to be searched if JBMAX is negative. 

For thermal stress problems(IEP>l) , FACLIM should be set to a 
number of the order one. 

Should be set to zero if JBMAX is positive; should be set to 
lower limit of desired range of eigenvalue to be searched if 
JBMAX is negative. 

SIGMA is not operative in the current version. Input of EPV and 
EPF on D card should now be used. 

For vibration problems set to some positive value. Vibration 
frequency will be eigenvalue times ANV. Leave blank for buckling 
problems. 

Print mode shape if MODE = 1. Leave blank if mode shapes not 
desired. Must be set = 1 to obtain plots. 

Prints total leads in all plates if LPLATE = 1. Leave blank 
if loads not desired. 

ALAMAX, RATLAM, EC, ANUC, DENSEC (4F10.4, F15.1, F7.3, F15.9) 

Panel length in x direction (See figure 1) 

Allows several solutions at different values of axial buckle 
length, A , to be. obtained. If ALAMIN is positive, A = 

ALAMIN x (RATLAM) 1 for i = 1,2,3. . until A greater than 
ALAMAX. If ALAMIN is negative A * EL/j for j = RATLAM, 

(RATLAM + 1) . . .ALAMAX 

If ALAMAX is negative a prescribed set of A's can be read on 
additional cards according to section 3C.4. of Reference 1. 

Young's modulus. Poison's ratio, and material density. If a 
number of plates are of the same isotropic material, these 
values will be used if none are input on the E cards. 

Applies only to the E cards. 


V 
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GROUP D 

IA, IB, IC, ID, IE, IG, I FAST , IEP , EPV , EPF, P (815, 3F10.6) 

IA Number of isotropic plates defined by group E cards. (There 

are 2*1 A group E cards). 

IB Number of orthotropic or anisotropic plates defined by group F 

cards (There are 3*IB group F cards). 

IC If IC * 0 no spring supports. If IC = 1 there are spring 

supports defined by a single group G card. 

ID Number of group H cards 

IE Number of group I cards 

IG Input of IG = 1 indicates material and layer properties will 

be read in. Otherwise IG = 0. 

IFAST Reduces computation time when calculating the buckling load at a 

number of wavelengths and only the minimum is desired. Should not 
be used for more than one eigenvalue at a given wavelength. 

0 no effect on calculations. 

1 Converged solution is determined only if lower than that 
for previous A. The number of negative roots at FACTOR = 

0 is checked. 

2 Same as IFAST = 1 except number of negative roots at 
FACTOR = 0 is not checked. One iteration is saved for 

each wavelength but if negative roots do occur at FACTOR = 0, 
solution may not converge. 

IEP Controls strain or temperature input. 

0 No strain or temperature input. 

1 Longitudinal loads calculated using strains EPV and EPF. 
Temperature may be present. 

2 Fixed temperature input - axial load to cause buckling 
is determined. 

3 Fixed load input P. Temperature distribution to cause 
buckling is determined. 

No E or F cards may be used for IEP = 2 or 3. Group DD cards can always be used. 
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EPV The variable longitudinal load NLV is calculated based on the 

strain EPV for any plate which NLV is read in as zero if 
IEP = 1. Used only when IEP ■ 1. 

EPF Same as EPV except it applies to fixed loading system. 

P Used only when IEP = 3. Value ot total fixed axial load on 

the panel . 

GROUP DP Input when IG > 0. These cards are mandatory when IEP = 2 or 3 
(no E or F cards). They are divided into 4 different categories: (l)Material 
properties, (2) Individual layer dimensions and temperatures, (3) Layer 
stacking sequence, and (4) plate definition and loading. 

DD1 Cards (Orthotropic Material Properties) 


E]_, E2 , 

E12, NU1 , RHO, ALFA1 , ALFA2 (7E10.3) 



El 

Young's modulus in 1 direction 



E2 

Young's modulus in 2 direction 



E12 

Inplane shear modulus 



NU1 

Poisson's ratio such that NU1*E2 

= nu; 

2*E1 

RHO 

Mass density of material 



ALFA! 

Coefficient of thermal expansion 

in 1 

di recti on 

ALFA2 

Coefficient of thermal expansion 

in 2 

di recti on 


Materials are internally numbered in sequence starting at 1 for the first card, 
2 for the second, etc. A card with a negative number in the first field 
signals that the previous card was the end of the DD1 cards. 

DD2 Cards (Orthotropic Layer Characteristics! 

MATL , H, ANGLE , T1V , T2V, T1F, T2F (15, 6E10.3) 


MATL Material number 

H Layer thickness 

ANGLE Angle between material axis and plate axis. See figure 1 

for positive direction. Must be in degrees. 

T1V, T2V Variable temperatures at inner and outer edges of layer, 
respectively. Inner eJge has least value of z. Variable 
temperatures are mu 1 ti plied by the eigenvalue. 

TIP , T2F Same as T1V, T2V except applies to fixed temperatures. 
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Layers are numbered in sequence starting at 1 for first DD2 card, 2 for second, 
etc. A card with a negative integer in the first field signals that the 
previous card was the end of the DD2 cards. 

DD3 Cards (Wall Generation) 

LAYNO(I ) (2014) I = 1 , 20 

LAYNO (I) Layer number of the Ith layer 

Layers are stacked in sequence in the direction of increasing z. The 
number of integers on the card will be equal to the number of layers for that 
particular wall. If the wall section contains more than 20 layers, continue 
on the next card by placing a zero in the first field followed by the 
additional layer numbers. More than one continuation card may be used so 
there is no limit to the number of layers in one wall. 

The first DD3 card is wall number 1, the second (if not a continuation 
card) wall number 2, etc. A card with a negative number in the first field 
signals that the previous card was the end of the DD3 cards. 

DD4 Cards (Plate Generation) 

JPA, I WALL , B, ANTV, ANLV, ANSV, ANTF, ANLF, ANSF (215, F10.5, 6F10.2) 

JPA Plate number such that 1 < JPA ^.120 

IWALL Wall number defined by DD3 cards. If it is desired to neglect 

anisotropic stiffness terms set IWALL to the negative of the 
wall number. 

B Plate width 

ANTV Inplane transverse, longitudinal, and shear plate loadings that 

ANLV are multiplied by the eigenvalue. If ANLV is input as zero and 

ANSV IEP = 1, ANLV is calculated from equation A7 using EPV. 

ANTF Inplane transverse, longitudinal, and shear plate loading that 

ANLF make up the fixed load system. If ANLF is input as zero and 

ANSF IEP =1, ANLF is calculated from equation A7 using EPF. 

A negative number in the first field signals that the previous card was the 
end of the DD4 cards. 

GROUP E (Isotropic Plates) 


Two consecutive cards, El, and E2 required for each plate. 
JPA, T, B, E, ANU, DENSE (15, 2F10.5, F15.1 F5.3, F15.9) 
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JPA Plate number 

T Plate thickness 

B Plate width 

E.ANU, Young's modulus, Poisson's ratio. and density for plate 

DENSE material. If any one is left blank value given on group C 

card will be used. 

E2 Cards 

ANTV, ANLV, ANSV, ANTF, ANLF, ANSF (6F10.2) 

ANTV Inplane transverse, longitudinal, and shear plate loadings that 

ANLV are multiplied by the eigenvalue. If ANLV is input as zero and 

ANSV IEP =1, ANLV is calculated from equation A7 using EPV. 

ANTF Inplane transverse, longitudinal, and shear plate loadings 

ANLF that make up the fixed load system. If ANLF is input as zero 

ANSF and IEP = 1, ANLF is calculated from equation A7 using EPF. 

GROUP F (Orthotropic and Anisotropic Plates) 

Three consecutive cards, FI, F2, and F3, reouired for each plate. 

FI Cards 

JPA, B, ANTV, ANLV, ANSV, ANTF, ANLF, ANSF (15, F10.5, 6F10.2) 

JPA Plate number 

B Plate width 

ANTV Inplane transverse, longitudinal, and shear plate loadings that 

ANLV are multiplied by the eigenvalue. If ANLV is input as zero and 

ANSV IEP = 1, ANLV is calculated from equation A7 using EPV. 

ANTF Inplane transverse, longitudinal, and shear plate loading that 

ANLF make up the fixed load system. If ANLF is input as zero and 

ANSF IEP = 1, ANLF is calculated from equation A7 using EPF. 

F2 Cards 

Dll, D12, 022, D33, D13, D23 (6E13.6) 

Dij Bending stiffnesses defined by equation A2 
F3 Cards 


AM, All, A12, A22 , A33 (F15.8, 4E15.8) 


jiiflttfe ! i ! ' ■ - : ; WSiKI !ti ■! . i i -i:i iu im .. : ...•b'iniMiMjitt: ! H li l 1 -■ i • i f ; V jjjj 
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AM Mass per unit area 

Ai.i Inplane stiffnesses defined by equation A1 

GROUP G (Spring Stiffnesses) 

Required (8E10.3) only if IC * 1. A single card which gives up to 8 
different values of spring stiffness that can be placed anywhere in the 
structure by boundary condition cards in Group H and I. 

GROUP H (Generation of Substructures) 

(2014) New plates or substructures can be defined from those 
previously generated by (1) rotation of a previously defined plate or 
substructure, (2) displacement of end ooints of previously defined plate 
or substructure, and (3) stringing together previously defined plates or 
substructures. 

Denote the number on Group H cards by i^, i 2 , i 3 . .. Than 
ij Number of additional numbers on the card. 

1 2 Number of plate or substructure being generated. 

1 3 Number of previously defined plate or substructure that is 
being modified. To accomplish a rotation enter i' 2 negative. 

Rotation .- To accomplish a rotation enter i 2 negative and i 3 positive. The 
counterclockwise rotation is i 4 x lt^S. (In this case, ij must equal 4.) 

Displacement .- To displace end points of a plate or substructure from node 
points, enter both i 2 and i 3 as negative. Then 

1*5 

EY1 = 1 4 x 10 3 

EZl=i 6 x ,’ 7 

ig 

EY2 = ig x 10 y 

1 1 

EZ2 « i 1Q x 10 1 

For this case i, must equal 10. EYj and EZj are the coordinates of the node 
points assuming the origin of the Y, Z coordinatj system is at the inter- 
section of the plate edge with its neutral surface. 


r 
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Multiplate Substructures and Boundary Conditions. - 

To string together a series of plates and substructure both ig and 
i 3 are entered as positive. Only substructures having nodes that 
form a simple chair, can be constructed. If a substructure has one end that 
is not attached to other structure it is a singly connected substructure and 
must be numbered starting at 121. The unattached end must be the initial 
point of this substructure. All substructures have nodes temporarily 
numbered 1,2. . . N y . Thus i 3 is the first plate or doubly connected 

substructure number of this sequence. If other plates or doubly connected 
substructures connect nodes 1 and 2: they are entered but with a negative 
sign. If singly connected substructures are attached at node 1 they are 
entered with a neaative siqn. Boundary conditions at node 1 are also entered 

for node 1 using a negative sign as -9xxxx00 where the x's indicate restraint 
against rotation. w, v, u, respectively. An x 3 0 means no restraint, x = 1-8 
selects the corresponding spring constant read on the G card and x = 9 is an 

infinite restraint. Specification of a free edge anywhere in the structure is 
not necessary. Such a specification will cause a program error. After all 
elements connecting 1 and 2 and boundary conditions and singly connected 
substructures at 1 are accounted for the next number will be positive 
indicating a plate or substructure connecting 2 and 3/again additional 
nega ive numbers indicate other connections to 3 or boundary conditions at 2. 
This process is repeated until the substructure is described. If it cannot 
fit on one card an intermediate substructure can be defined and used to start 
a new card. 

If a singly connected substructure or a boundary condition is at the 
last node (N ) enter a zero then the substructure or boundary condition 
preceded by a negative sign. 

For singly connected substructures, N must be greater than one, and for 
doubly connected substructures N must be yreater than two. Violation of 
this rule will cause i neon act results or a dump. A diagnostic in the 
plotting program will be printed if this rule is violated. 

Examples of the generation of a variety of substructures are illustrated 
in the sample problems of reference 1. 

GROUP I (Generation of Final Structure) 

(2014) 

Final structure is assembled from plates and substructures in Group I 
cards with format 14. If as in Group H we denote numbers on Group I as 
ip i* 2 - . then 

i-j Number of additional numbers in the card. 

i 2 Node number. One card must be prepared for each node connected 

to a higher numbered node or having a prescribed boundary 
condition. 

i *3 Node number (must be greater than i‘ 2 ) that is connected to i 2 

i^ Plate or substructure that connects the two nodes. 
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1r“i Additional pairs of numbers, the first being the number of 

1 the node which is connected to i 2 and the second, the 

number of the plate or substructure making the connection. 

Note, if more than one element connects the two nodes the 
higher numbered node number will be repeated. 

i n+ 1 Boundary condition or the number of a singly connected 

substructure at node i 2 preceded by a minus sign. Notation 

for boundary conditions is the same as in the H cards 
(9xxxx00). 

i n+9 Additional singly connected substructure or a boundary condition 

but not preceded by a minus sign. 


EXAMPLE PROBLEMS 


The first problem is buckling under compression and shear of the panel in 
figure 3. it is desired to obtain the two lowest eigenvalues for values of 
buckle length A from 2 to 32. Buckling modes, including plots, are desired. 

The input data using the option for layered material property input, .is shown 
ii Table 1. The shear flow in the hat and skin was calculated by equating 
the shear deformation of the two paths accounting for the inplane shear 
stiffness. A portion of the output is shown in Table 2. Plots of the 
structure, buckling mode shape, and a plot of buckling strain versus X are 
shown in figures 4 to 6. Note in figure 4 node locations are shown by 
circular symbols. 

Several capacities of the program are illustrated in this example. 

The flange of the hat section is bonded to the skin so that the two act 
together. Such a plate results in a non zero B matrix which cannot be 
handled in VIPASA. However, it still can be treated accurately by the 
procedure given in reference 3. The plate is modeled as a substructure with 
both the skin and flange tied together at a spacing of 1/8 the total plate 
width. The flange is given an eccentricity of one half the sum of the two 
plate thicknesses to give the resulting substructure the proper bending 
stiffness. 

Another feature is the capability to treat structures with periodic 
buckling patterns in the transverse unction by giving the boundaries the same 
node number. Note in figure 4 that both ,'dges of the panel are numbered 5 which 
causes the deformation on the left edge ti be compatible with the deformations on 
the right edge. At the end of table 2, the •'iagnostic is printed that node 5 is 
in error. This is the result of assigning two different locations to the same ncde. 

In order to construct the plate going from node 1 to node 5 it is 
necessary to rotate a plate by 180°. Even though this plate is physically 
identical to plate number 5, plate 5 cannot be used since it would not have 
the shear force in the proper direction or the fibers in the proper sequence 
upon rotation. Therefore, a new plate is generated with tr e sign on the 
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angle plys reversed and the shear load reversed so that it will be identical 
to plate 5 after a 180° rotation. 

The thermal stress capability is illustrated in the second problem. 

It is assumed the panel of problem I was in a stress free state at an elevated 
curing temperature. The effect of a uniform temperature drop of 300° F is 
investigated with option IEP = 2. The input cards are shown In table 3 with 
portions of the output in table 4. The thermal stress in each plate number 
is shown for no axial load and the prestress assumed for the buckling analysis 
is shown. The the ;nal stress was rather small and. had little effect on the 
buckling load. 

A check of the accuracy of the program for thermal stress problems is 
indicated in the third example. The compressive buckling load of a nonuniform 
temperature, titanium zee stiffened panel which was determined in reference 4 
has also been determined with VIPASA. The input for the problem is given in 
cable 5 and the pertinent output in table 6. Results from reference 4 are also 
shown for comparison. There is generally good agreement; the small difference 
in buckling load is attributed to effects of eccentric connections which were 
ignored in VIPASA. 


COLLUDING REMARKS 


Extensions developed at Langley Re^?arch Center to the VIPASA computer 
program for panel buckling and vibration analysis have been described. A 
condensed version of the users manual including new input options is given. 
Example problems with sufficient documentation for checking results are also 
given. 
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APPENDIX A 


STRESS STRAIN RELATIONSHIPS AND THERMAL STRESS CALCULATION 

The relations defining the elastic properties including temperature 
effects are: 
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where 



(A3) 


C jT * f Q n Tdz (A4) 

See figures 1 and 2 for positive values of variables. For an orthotropic 
lamina inclined at an angle® from the plate coordinate system the Q's are 
first calculated in the material coordinate system and then transformed to 
the plate coordinate system as follows: 


Q' 
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(A5) 
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^11 = n cos4e + 2(Q' 12 + 2Q' 33 )sin 2 6cos 2 9 + Q' 2 2 sin4e 
Qi2 = ^'n + Q '22 ’ 4Q^3 )sin 2 0cos^e + Q'^fsin 4 © + cos 4 0) 

Q 22 = Q'iiSin 4 0 + 2 (Q * 12 + 2Q' 33 )sin 2 0cos 2 0 + Q' 22 cos 4 0 

Q 13 = ^±11 “ Q 'l2 “ ?Q'3 3 )sin0cos 3 0 + (Q' 12 - Q' 22 + 2Q' 33 )sin 3 0cos0 

Q23 = (Q'u “ Q' 12 " 2Q' 33 )sin 3 0cos0 + {Q ' 12 - Q' 22 + 2Q' 33 )sinecos 3 0 (A6) 

Q33 = (Q' n + Q' 22 - 2Q' 12 -2Q' 33 )sin 2 0cos 2 e + Q' 33 (sin 4 0 + cos 4 0) 

QjT = Q' 1T cos 2 0 + Q' 2T sin 2 e 
Q2t = Q'ijSin 2 9 + Q' 2 t cos 2 6 
Q 3T = (Q' 1T - Q' 2T )sin0cose 


Positive values of 6 are in the direction of rotating x into y. 


The limitations of viPASA require that A, 3 = Ag, = 0 as well as the 
entire B matrix must be zero. Symmetric balanced laminates will always 
satisfy this condition. When layer properties are input the A, B, and D 
matrices are printed and should be examined to insure input is valid. 
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With the above assumptions the transverse strain may be eliminated 
from the equations to yield 

N L * (A 11 * A 12 2 /A2 2 )(-?u/3x) + A 12 (N t - C 2T )/A 22 t C 1T (A7) 


The total load on the panel cross section Is 

P = -au/SxJj^b. + X>1 b l 
The strain required to achieve this load Is simply 


-wax =■ - (A9) 

2>i b i 

where 

a i = A 11 " A 12 '^2 for plate 1 

b.. width of plate i (A10) 

c i = A 12 (N T‘ C 2T )/A 22 + C 2T f ° r plate 1 

Equation A7 is used to calculate the longitudinal loading on all the elements 
for a given longitudinal strain. For thermal stress calculations two options 
are possible. For a fixed temperature input (IEP = 2) the axial load required 
for buckling is determined by assuming the fixed longitudinal strain, EPF = 
-.02 and the variable strain EPV = .01 . (The large fixed tensile strain is 
to insure that buckling does not occur under the fixed loading conditions). 
Prior to this calculation the prestress is evaluated for a strain determined 
from equation A9 with P = 0. This solution is printed first followed by the 
prestress solution used for the buckling analysis. 

For a variable temperature input with a fixed axial load P (IEP = 3). 

The fixed strain is calculated from equation (A9) and the variable strain also 
from equation A9 with P = 0. Thus the fixed load system results in an axial 
load P while the variable system results in no axial load. The eigenvalue 
multiplied by the variable temperature input then represents the temperature 
distribution required for buckling i,. the presence of the fixed axial load. 
Again the calculated prestress solution is printed. 


/ 
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It can be seen that the assumptions of constant strain over panel cross 
section is not good in all cases so caution should be exercised in the 
application of the thermal stress analysis. However, in some cases such a 
simple analysis is adequate, such as symmetrically heated cross sections or 
panels over many supports so that thermal bowing is restrained. 
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APPENDIX B 


LANGUAGE: 

PURPOSE: 


USE: 


DESCRIPTION OF PLOTTING SUBROUTINES 


SUBROUTINE PSEUDO 


COMPASS 

Tc create and write an appropriately named Plot Vector 
File. Through linkages set up hy an initial call to 
PSEUDO, all subsequent graphics data generated by the 
user will be routed through one of the PSEUDO entry 
points and written on the Plot Vector File. The PSEUDO 
processor is designed for use with the frame dependent post- 
processors described in Section 1-3* Volume IV, of the 
Computer Programing Manual. 

CALL PSEUDO 

or 

CALL PSEUDO (FN) 

FN tile name left- Justified with zero fill. 

Default file name is SAVPLT. 


Example : 

CALL PSEUDO 

This will establish a Plot Vector File named 
SAVPLT. 

CALL PSEUDO (6LMYFILE) 

This will establish a Plot Vector File named 
MYFILE. 

NOTE : The Plot Vector File (or Files) will usually be 

written to disk (as opposed to tape) and may be 
postprocessed following user program termination 
via appropriate specification of one or more 
PLOT control cards (see Section 1.3, Volume 
IV, Computer Programing Manual). 

(l) An initi.ilizing call to PSEUDO (with or without a 

file name argument) must be made prior to any calls 
to CALPLT or any other graphics output routine. 


RESTRICTIONS: 


COMPUTER 
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PSEUDO 


(2) Every Plot Vector File should be terminated 
with a 999 pen code, CALL CALPLT(0. 0,0. 0,999). 
The transmission of the 999 code will cause an 
EOF write on the Plot Vector File, and the file 
will temporarily be closed. Thus, any given 
Plot Vector File will contain only one 999 pen 
code and/or one EOF. 

(3) To continue plotting execution following trans- 
mission of a 999 code to a current Plot Vector 
File, the user program must call the PSEUDO pro- 
cessor to create new Plot Vector File (i.e., 

CALL PSEUDO(6LMYFIL2)). 


METHOD: In addition to entry PSEUDO, this processor contains two 

other entry points, namely PLT9999 and PLT9998. An initial-| 
izing call tr PSEUDO will set FLT9999 into the processor 
switching mechanism (PLOTSW). Subsequent plot data gen- 
eration will then be routed via CALPLT, PLOTSW, and 
PLT9999 and written on tie Plot Vector File. The entry 
PLT9998 is used to record special purpose data from 
routines NFRAME and PLTS f1 0P. 


ACCURACY: 

REFERENCES: See Section 1.3, Volume IV, Computer Programing Manual. 

STORAGE: 2155q locations total for direct subprograms 

SUBPROGRAMS USED: NUMARG, PLOTSW 

OTHER CODING 
INFORMATION: 

SOURCE: E. C. Johnson, NASA-Langley Research Center 

QUESTIONS ON THE USE OF THIS PROGRAM SHOULD BE DIRECTED TO THE ACD PROGRAMER 
SUPPORT GROUP, EXTENSION 35^8. 
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SUBROUTINE NFRAME 


FORTRAN 

To provide users specific means of executing frame 
advance movements on any plotter device via an appro- 
priate frame oriented device postprocessor. Frame 
advance distances are generally defined to be incre- 
mental from current frame origin (i.e., comparable to 
frame advance executions for the DDI or 252 CRT devices) . 
CALL NFRAME is intended to be used as a frame advance 
mechanism, n£$ as a plot origin offset. 

CALL NFRAME 


CALL NFRAME(H,V) 


where 


H and V 


are Horizontal (parallel to device X) and 
Vertical (parallel to device Y) distances 
from the current frame origin. H and/or 
V must be expressed in floating point 
inches. 


The short form CALL NFRAME will cause the device postpro- 
cessor to execute a frame advance move parallel to the 
device X (horizontal) axis. The movement will be the 
maximum horizontal distance used in inches plus h inches, 
where h will be an increment appropriate to the particular 
device (CXh*^"). The frame advance will be rounded to 
whole inches. 
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VOLUME f V 
SECTION 1.4.3 
ORG. 11.910 
DATE \/ 1/75 


COMPUTER 


M AN UAL 


COPY DATED none 


PAOC 1 or 7 


SUBROUTINE CALPLT 

LANGUAGE: FORTRAN 

PURPOSE: To move the plotter pea to a new location with pen up 

or down. 


USE: 


CALL CALPLT (X,Y,IPEN) 
where 


X,Y 


are the floating point values for pen move- 
ment. 


IPEN = 2 pen down 


= 3 pen up 

Negative IPEN will assign X = 0, Y = 0 as 
the location of the pen after moving the 
X,Y (create a new' reference point). 

= 999 Writes a terminating block address of 999 
to terminate the Plot Vector File and all 
further processing is skipped. 

CALL CALPLT(0.0,n. 0,999) 


RESTRICTIONS: All X and Y coordinates must be expressed as floating 

point inches (actual page dimensions) in deflection from 
the origin. 


A TERMINATING BLOCK ADDRESS (IPEN = 999) MUST BE GIVEN 
AS THE LAST PLOTTING INSTRUCTION BEFORE ENDING A PROGRAM 
WHICH USES ANY OF THE PLOTTER SUBROUTINES: THIS IS TO 

BE SURE THAT ALL PLOTTER INSTRUCTIONS ARE WRITTEN ON THE 
PLOTTER TAPE. 
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VOLUME IV 
SECTION 1 *4.35 

COMPUYEi 
PROP R AMINA 
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SUBROUTINE NOTATE 


LANGUAGE : FORTRAN 

PURPOSE: To draw alphanumeric information for annotation and label- 

ing and provide special centered symbols for annotation of 
data points. 

USE: CALL NOTATEfX.Y.HEIGHT^CD.THETA.NOCHAR) 

where 

X,Y are the floating point coordinates of the 

first character. 

For alphanumeric characters, the coordinates 
of the lower left-hand corner of the characters 
are specified. 

For special symbols 0 - 5» the coordinates of 
the center of the symbol are specified. 

For special symbols above 6, the coordinates 
of the lower left-hand corner of the character 
are specified. 

HEIGHT specifies character size and spacing in float- 

ing point inches for a full-size character. 

The smallest possible character is 0.07 inch 
high. The width of a character will be (L/7)* 
HEIGHT and the space between characters is 
(2/7 )*HEIGHT + WADJ (see Figure l), where 
WADJ is width adjustment. 

The ith character is plotted at: 

Xi = X + (i-l)( 6/7) (HEIGHT) (COS 0) + 

(i-1) (WADJ) (COS 9) 

y ± * Y ♦ (i-1) (6/7) (HEIGHT) (SIN 9) + 

(i-1) (WADJ) (SIN 0) 

where WADJ = 0 for HEIGHT > .1 

= .01 for HEIGHT <l .1 
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H — H 


H a HEIGHT 

W = (U/7) « HEIGHT 

S = (6/7) * HEIGHT + WADJ 


Figure 1 


i s the string of characters to be dT*awn and is 

usually written in the form: nHXXXX (the 

same way an alpha message is written using 
FORTRAN format statements). Instead of speci- 
fying alpha information as above, one may give 
the beginning storage location of an array con- 
taining alphanumeric information. 

Special symbols will be drawn when BCD is an 
integer i ference and NOCHAR it. negative (see 
Figure _,/. 

NOTE : A binary zero in the BCD string will 

cause truncation of plotting at that 
point , and a normal return to the call- 
ing program. 
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NOTATE 


h 


THETA is the angle in floating point degrees at 

which the information is to be drawn. Zero 
degrees will print horizontally reading from 
left to right, 90° will print the line verti- 
cally reading from bottom to top, 180° will 
print the line horizontally reading from right 
to left (i.e., upside down), and 270 ° will 
print vertically reading from top to bottom. 

NOCHAR is the number of characters, including blanks, 

in the label. A negative NOCHAR will produce 
a single special symbol from the integer ref- 
erence table. (See METHOD for further explana- 
tion. ) 

RESTRICTIONS: Noted under METHOD. 

METHOD: The character height is a variable entry parameter to the 

subroutine NOTATE. However, the width-to-height ratio is 
fixed at 4/7. This is because the characters are defined 
by a series of bi-octal offset pairs for a by 7 matrix 
as showr by the examples in Figure 2. The reference origin 
for the offset pairs which define each character is the lower 
left-hand comer of the matrix. The X and Y values 
which are entry parameters to NOTATE define the location of 
the lower left-hand corner of the first character to be 
plotted ■'or this entry to NOTATE. Subsequent characters 
to be plotted are spaced from the previous character origin 
by 6/7 of the specified character height plus width adjust- 
ment. 
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SUBROUTINE PNTPLT 


LANGUAGE: FORTRAN 

PURPOSE: To draw NASA Standard Plot symbols centered on a given 

coordinate value. 

USE : CALL PNTPLT (X, Y, ISYM, IS) 

where 

X is the X coordinate for the centered symbol 

in floating point inches. 

Y is the Y coordinate for the centered symbol 

in floating point inches. 

ISYM is an integer specifying the symbol to be used. 
(See Figures 1 and 2.) 

- 21 for a point * 

= 22 for a plus sign + 

IS is an integer value specifying the size symbol 

to be used. 

= 1 small 

* 2 medium 

=* 3 large 

(See Figure 1.) 

RESTRICTIONS: 

METHOD: 

ACCURACY: A positive integer value for ISYM ju the calling sequence 

will produce symbols of the same quality as in Figure 10. 
A negative integer value will produce symbols of less 
quality but will result in a considerably faster computer 
run. 
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NASA STANDARD PLOT SYMBOLS 

INTEGER SIZE 

REFERENCE SMALL MEDIUM LARGE 


11 

® 

© 

© 

12 

IS 

m 

E 

13 

<3> 

<e> 

<$> 

14 

A 

A 

A 

15 

fk 

k 

k 

16 

& 

& 

& 

17 

© 

© 

ffl 

18 



<5> 

19 

♦ 

♦ 

❖ 

20 

& 

rib 

A 


FNTPUT 

i 


Figure 2 
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SUBROUTINE NUMBER 


LANGUAGE: 

PURPOSE: 


USE: 


RESTRICTIONS: 

METHOD: 

ACCURACY: 

REFERENCES: 


FORTRAN 

To convert a floating point number to BCD (expressed in 
F format), and draw the resulting alphanumeric charac- 
ters. 

CALL NUMBER(X,Y,HEIGHT,FFN,THETA, NODIGIT) 
where 

X,Y are the coordinates in floating point inches 

of the left lower corner of the first digit 
of output. 

HEIGHT is the height of the plotted number in float- 
ing point inches (see NOTATE routine). 

FPN is the floating point number to be drawn. 

THETA is the angle in floating point degrees at which 

the number is to be drawn (see NOTATE routine). 

NODIGIT is the number of decimal digits to the right 
of the decimal point for output. 

N0DIGIT=-1 or N0DIGIT=0 both specify no decimal 
places; however, -1 suppresses the decimal point. 

The number is restricted to a maximum of 12 significant 
digits. 


The routine trunc&ies the floating point number at the 
required decimal place. 




NASA 


LANGLEY 

RESEARCH 

CENTER 


VOLUME 1 V 
SECTION 1.4,40 
ORO. 11.510 
DATE 1 / 1/75 


REPLACES COPY DATED NONE 



SUBROUTINE AXES 


LANGUAGE: 


PURPOSE: 


FORTRAN 

To draw a line, annotate the value of the variable at 
specified intervals with or without tic marks, and pro- 
vide an axis identification label. 

CALL AXES (X,Y, THETA, S,ORG,SFX,TMAJ,TMIN, BCD, HEIGHT, NOCHAR) 


where 


THETA 


are the coordinates in floating point inches 
of the starting point of the axis with refer- 
ence to the plotting area origin as established 
by CALPLT . 

is the angle of rotation measured counter- 
clockwise from the X-axis in floating point 
degrees. Normally, THETA is 0° for an X-axis 
and 90° for a Y-axis. 

is the length of the axis in floating point 
inches. Should be a multiple of TMAJ. 

+S will generate tic marks. 

-S will eliminate tic marks. 

is the functional value to be assigned to the 
origin (i.e., the value of the first scale) in 
floating point. 

is the adjusted scale factor for the array to 
be plotted ( change in value per inch ) . 

NOTE: Values of ORG and SFX which will produce 

a reasonable scale may be calculated using the 
subroutine AS C ALE or BSC ALE. 

is the distance in floating point inches for 
major tic marks (0.25 inches high). Numbers 
are placed on the axis at the major tic marks 
in accordance with the values of ORG and SFX. 
The numbers written along the axis are adjusted 
to be between 1000.00 and 0.01 in magnitude. 
Immediately after the last number on the axis 
is placed the caption xlO^xp, where exp is 
the required exponent. 
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If the values are integer multiples, the deci- 
mal point and decimal places are eliminated. 

A negative TMAJ will cause the actual value to 
be written instead of the adjusted value. 

TMIN is the number of divisions per inch in 

floating point for minor tic marks (0.12!? 
inches high). To eliminate minor tic marks 
the following may be used: 

TMIN = 0. 

BCD is the character label for the axis (see 

NOTATE routine). 

HEIGHT is the height of the full-size characters in 
the BCD title. Numbers at the tic marks will 
be (0.75 -^HEIGHT) '.r^h. HEIGHT is in floating 
point inches. 

If HEIGHT — 0. , all annotation will be eliminated. 

NOCHAR is an integer specifying the number of charac- 
ters in BCD title. A negative NOCHAR places the 
annotation on the clockwise side of the axis 
and a positive NOCHAR places the annotation on 
the counterclockwise side of the axis. NOCHAR = 0 
is not allowed. If it is desired to have no 
label, then the BCD parameter should be 1H , 
and NOCHAR either +1 or -1. 

RESTRICTIONS: Only perpendicular axes are recommended. 

METHOD: 

ACCURACY: 


REFERENCES: 





Table 1.- Input for Example Problem Number One. 


28 ?b 13 2QJ) 4J0 . 2P _ . 2 


PPOBLEf 4 3.C UD • S. F LAfifi ESI . S.. H .£AKg <t5.iL> 

2 1 1QQQQ..Q Z* 



L-11E6.. 


32.1.4l421354i 

4 1 0 1 

.4ftBEt> _ *19. 


.•iUU. 

• 05.7 ... 


_fL» 0 ... ILail 


. J 052 0.*_Q. 

1 «JLQJ52 


1 .0052 

45.0 







1 .0104 

(LlQ ... 







1 .052jL._ 

0.0.. 







-1 









2 3 1 

... 2 3 

2 

3 

4 

. 3 2 

2.. 

2 l 

.3 2 

P 3 1 
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__2 ... 

2 

__L 

3_ 2. 




2 2 .... 1_ 
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_2 

... .1 3. 
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2 _2 1 . 

? 3. 

.2 _ 

3 

_3.. 
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1. 3 

2 . . 

3 2 L. 

• 3 2. 

. 2. 

2 

_iL 

3 . 2 

_ 2._ 
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3 

- 1 









1 4 

0.125 








2 4 

2.580 








3 2 

1.32b 








4 3 

1.250 








5 4 

1.310 








6 5 

1.310 








7 2 

0.125 








-1 









4-9 3 

300 








4 -10 _ 3_ 

80 








4 1L 9 

4 10 








3 12 5 

5 








10 -13 -7 
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-4 
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0 -624 

-4 



4 tJL4__. 6 
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5 8 1 

-13 1 
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5 3 8 

8 8 

8 







5 8 3 

2 -11 

3 







5 1 2 
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14 . 
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1? 
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334*84_ 

264,27 

I85.J3 

185.73 

450*00 

-450*00 
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Table 2.- Selected Output from Example Problem Number One, 

VIPASA RESULTS _ _ 

INPUT DEVICE NUMBER Ts 5 OUTP U T~DE V ICE NUMBER IS ~6 
NMK NSR NSC I PAST NSN NSK NOLAM 

28 26 13 20 0 400 2 20 


I PLOT NSPAC PAGE 

2 -0 - 0.00 

NUMBER OF PROBL EMS IS 1 

DATA FOR PROBLEM NUMBER - 1 




PROBLEM 3C 
2 1 
32.0000 


(0_*S, F LAN GES) SHEARs450 . 

10000,0 2.000000 

2,0000 32.0000 1.414? 



0.000000 0.0000 0.0000 1 
- 0.0 - 0.000 - 0,00000 000 0 


the ABOVE TWO DATA CARDS HAVE THE FOLLOWING MEANING 

THE PROBLEM IS TO FIND ALL EIGENVALUES FROM NUMBER 1 DOWN TO NUMBER* 

INCLUSIVE. TO AN ACCURACY OF AT LEAST ONE IN 10000^0 

THE LFNGTH OF THF STRUCTURE IS ’ 32.000000 AND LAMBDA IS INCREASED FROM 
2.000000 TO THE HIGHEST POSSIBLE VALUE OF LAMBDA NOT EXCEEDING 
32.000000 BY SCALING SUCCESSIVELY BY 1.414214 

MODFS ARE FOUND AT ALL EIGENVALUES 

THE NUMBERS OF PLATE'S DEFINED IN SIMPLIFIED WAY, OF PLATES DEFINED _IN " 
GENERAL WAY. INDICATOR OF GENUINE SPRING SUPPORTS. NUMBERS OF 
SUB-STRUCf URES AND OF C A Rf) S_DE FINING FINAL STRUCTURE' 

0 0 0 9 4 


ARE, RESPECTIVELY, 

.. 


IG= 1 JEAs>T= 0 IEP= 1 

THE FIXED LONGITUDINAL STRAIN IS 
THE VAFLAfikEJ-ONGJ TUJHNAJL STRAIN 


- 0.000000 
IS .001000 


— 


MATERIAL INPUT 

"V 


MA.TLRIAL___E1 £2. El 2 Nyj 

1 . 1 9SE ♦ 08 . 1 1 IE *07 . 488E ♦ OB ,190E*00 .570E-01-0. 



. K lc;, . 


,,... . j i ; .,-• , ; • « * - tWr?.«V 



Table 2.- Continued. 


Q Q 

^ pj 
>d & 

8g 

w j5 



W ALL H ATL H A NGLE JT1V _T2V T 1 F T2F 



1 1 

.5200E-02 .4500E+02— 0 • 

-o. 

-0. 

• 

o 

* 

1 

•5200E-02 — • 4500E* 02— 0 • 

-0. 

-o.. 

-0. 

1 

•5200E-02 0. -0. 

-0. 

-o. 

-0. 

1 

.5200E-02 .450QE+02-0. 

-0. 

-0. 

-0. 

1 

.5200E-02 -.4500E+02-0. 

-0. 

•o. 

-0. 

1 

•5200E-02 .4500E+Q2-0* 

-0. 

-0. 

-0. 

1 

.5200E-02 -.4500E* 02-0 • 

-0. 

-0. 

-0. 

1 

.1040E-01 0. -0. 

-0. 

-0. 

-0. 

1 

.5200E-02 -.4500E*Q2-0. 

-0. 

-0. 

-0. 

1 

• 5200E-02 • 4500E + 02— 0 • 

-o. 

-o. 

-0. 

1 

.5200E-02 -.4500E*02-0. 

-0. 

-0. 

-0. 

1 

.5200E-02 . 450 0E+ 02-0 • 

-o. 

-0. 

-0. 

1 

• 520 OE-02 0. -0. 

-0. 

-0. 

-0. 

1 

• 5200E-02 - .4500E* 02-0 • 

-0. 

-0. 

-0. 

1 

•5200E-02 .4500E+02-0. 

-0. 

-0. 

-0. 


ZREF FOR 

WALL NO. 1= .4160E-01 





M ASS PEP UNIT AREA OF WALL NO. 1= .4742E-02 


A-HAT RIX FOR WA LL no. 1 P-MATRIXFOR NAL L N O, 1 


• 7657E*06 

•3027E+06 

0. 

•3952E+03 

• 1900E* 03 .3368E *02 

•3027E+06 

• 3824E*06 

0 • 

.1900E4-03 

•2362E+03 .33686*02 

0. 

0. 

•3257E+06 

•3368E+02 

• 3368E+02 .2033E+03 


THERMAL LOADING TERMS 

B-MATPI X FOR WALL 

NO. 1 

VARIABLE 

FIXED 


til 64E-09 

-• 1746E-09 

-.1455E-10 

0* 

0. 

-.1746E-09 

-.2910E-09 

-•7276E-1 1 

0. 

0. 

- • 1 455E-1 0 

-•7276E-11 

— .17 46E— 09 

0. 

0. 




•*.* tin r»^: r 


r fxltitiir. t it-i* .■'*«?• 


■ . • . • 


* 

' h 

Table 2.- Continued. 



HHH 



' , . 



NEGATIVE. MALL NUMB£fiS^ _I«CLLCAT^_ANISQ.TROPIC ST.I/_rNtSS_J£8»iS-_NEGLECTED_. 












Table 2.- Concluded 



EIGENVALUE NO. LAMBDA FACTO* 

1 . 3 1 99999*# (>E ♦ 0 2 . 29836 1 9 1 BE ♦ 0 1 • l jui > >ui. ^ u i • *. » 

T HE TOTAL AXIA L LOAPfS . 31 3205 08E +05 ___ __ 

THE MODF IS 

CLOCKWISE ROT. DEFLN. DOWN DEFLN. TO RT. LONG. DEFLN. 

REAL PARTS ON ALTERNATE LINES. WITH IMAGINARY PARTS INSET BENEATH 
. 1 0000000E+00 -.25599137E+01 . 8287 1 725E-0 3 . 75670988E-0 1 

.00352492 -.01483189 .00109697 .00013147 

-.99999997E- 01 - . 25 599 1 37E ♦ 0 1 - . 828 7 1 7 35E -0 3 _ . 75670983E-0 1 _ 
.00352492 .01483188 .00109697 -.00013147 

.99999997E-01 -. 25599 1 37E ♦ 0 1 .82871 735E-03 . 75670988E-0 1 

.00352492 -.01483188 .00109697 .00013147 

-. 1 OOOOOOOE^OO - • 25599 13 7E ♦ 0 1 -.82871724E-03 . 75670988E-0 1 

.00352492 .01483189 .00109697 -.00013147 

-.1 57444 25E-10. - . 2644 3 1 1 0E »0 1 . 1 6 1 0359 1 E- 1 1 .74481098E-0JI 

-.01996355 .00000000 .00110327 0.00000000 


LIMITS 


,298361918E*01 


ON FACTOR 


•298364258E*! 


.75670983E-01 

-.00013147 

.75670988E-01 

.00013147 

•75670988E-01 

-.00013147 

,74481098E-0_1 

0.00000000 



' 

J 




J . 

: 


ALL REQUIRED EIGENVALUES HAVE BEEN FOUND 

COORDINATES UF NODES_ 

NODE Y Z 


LAMBDA = 


32.000000 


1 _o* 

2 .4 5800 000F + 01 

3 _r 7200_0Q00E*01 

4 . 11780000 E *02 

5 . 13090000E*02 
Y COORDINATE OF NOUE= 5IS 






















Q *T 

Ei 


Table 4.- Selected Output from Example Problem Number Two. 


DATA 


PROBLEM NUMBER 


7 


PROBLEM ^3C 
1 1 
32.0000 


(O.S. FLANGES) SHEAB=ASO» TEMP=-30„. 


1 1 l^OOD-O 2.000000 OaOOOOOO 0.0000 O.OQM 1 

12.0000 2.0« 0 32* 0000 l.*»14Z -0.0 -0.000 -0.000000000 

ABOVE TWO DATA C_ARUS HAVE THE FOLLOWING MEANING 

THE PROBLEM IS TO FIND ALL EIGENVALUES FROM NUMBER 1 DOWN TO NUMBER 1 

INCLUSIV E , T O AN ACCURACY OF AT LEAST ONE JN 10000.0 

THE LENGTH OF THE STRUCTURE _ I S‘ 32.000000 AND LA 'BOA IS INCREASED FROM' 
2.000000 TO THE HIGhtST POSSIBLE VALUE OF LAMi DA NOT EXCEEDING 
32.000000 BY SCALING SUCCESSIVELY BY 1, 41421 4 _ " 

MODES ARE FOUND AT ALL EIGENVALUES 

THE NUMBtnS OF PLATES DEFINED IN SIMPLIFIED WAY, Or PLATES DEFINED IN 


0.0000 0.0000 1 0 

000 - 0.000000000 


DOWN TO NUMBER 



GENERAL WAY. INDICAT OR OF GENUINE S PRING S UPPORTS, NUMBERS OF 

SUB-STRUCTURES AND OF CARDS DErlNlNG F I NAL STRUCTURE ARE, RESPECT! 
0 0 0 9 4 _ 

_ _ 

= l IF AST = 2 IEP= 2 

BUCKLING LOAD WILL BE CALCULATED FOR F I Xt U TEMPERATUKE INPUT 

EIGENVALUE DETERMINED ONLY If_LESS THAN THAT FO R PRE VIOUS LAMBD A 
CONVERGENCE MAY NOT BE ObT A I NED IF THERE ARE NEGATIVE ROOTS AT F ACT OR = 


NUMBERS OF 

E ARE, RESPECTIVELY* 


MATFPIAL INPUT 

— . 








/■ 


I 



Table 4.- Continued. 


. 1 *5.?Q0E-0<L . 

1 . 52 QQ_t- 02 . - 

1 «5200E-_02 0. 
•5200E-02 


L .520 0E-02 - 

1_ .5200E-02 

Il_ .520 0E-0 2 - 
A _i104M.-0T.JQ 

„ L_. .5200E-02 - 

_1 «52 QQ g .rQ £. - 

I .5200E-02 - 

l_ *520 0E-02 
A _ •52pOE^Oj2 .0 
l„__*520 0E-02 - 
1 .5200 E- 02 


JMLL._ 

Jl 4500E+_02« 

* * 

------ 


t 4500E_*U2' 

*4500E*02« 

. 4500 E*02» 

* ' 

.« 4500 E+ 0 j 2 * 
■ 45 0 OE* 0 2 « 
a^50lOE ♦0 2* 
*4500t+02- 
• * 

• 4500fc4 02^ 
. 450 OE* 02« 


ZREF FOR WALL NO. 1= .MbOE-Ol 

MASS PER UNIT AREA OF WALL NO. ' t « 

A -MATRIX FOR WALL NO. 1 

, 7657 E* 06 • 3^027E + 06 0 . 

.3027E + 06 • 3824E* 06 0.„ __ 

0. 6. • 3257E + 06 


_ 8-MATRIX FO R WALL NO. I 

• 1 164E-09 - • 1746E-09 !L.14b5EM0 
- . 1 746E-09 -.2910E-09 -.7276E-H 
-.1455F-10 -.7276E-1 1 -.1746E-09 


T2V " T 1 F T2F 

- , 3000E *03 -• 3000E + 03 
-.3000E403 -.3000E*03 
-.3000E + 03L -.3000E.403 
„V*30qOE403_-.3000JE403 
“ ^.30001403 -.30b0E*b3 

- .3 00 OE ♦ 03_ - • 3 0 OJ) E ♦ 03 
- .3000 E4Q3 -.3 000 E403 

“>,3000E403 -#_30J>1E4 03 

- • 3 0 00. E ♦ 0 3__- ._3 0 0 0 E* 03 
-.30^0E*03 -.300b _ E403_ 
-.30005403 -.30002403 
-.3000E403 -.3000E403 
-.3000E403 -.3000E403 
-.3000E403 -.3000E403 
-.3000E403 -.3000E4_03_ 


•4742E-02 

D-MATHIX FOR WALL NO. 1 

. 39 52 EV03 '.190 0E> 03 ". 3 366E 4*0 2 
• 1 900E + 03 .2362E403 .3368E402 

.3368E402 .3368E402 .20331403 

THERMAL LOADING TERMS 
* VARIABLE F jfXED _ 

... -• 2977 E 493 

0. -.3458E403 

' b. 0 . 





Table 4.- Concluded 


THERMAL ST RESS F O R AX IAL L OAD P « 0* 


PLATE 

number breadth 

ntv 

NLV 

NSV 

NTF 

NLF 

NSF 


1 

•12500 

-0,00 

0.00 

-0.00 

-0.00 

-12.24 

334.84 


2 

2.58000 

-0.00 

0.00 

-0.00 

-0.00 

-12.24 

264.27 


3 

1,32800 

-0.00 

0.00 

-0.00 

-0.00 

-1.29 

185.73 


4 

1.25000 

-0.00 

0.00 

-0.00 

-0.00 

75.30 

185.73 


5 

1,31000 

-0.00 

0.00 

-0.00 

-0.00 

-12.24 

450.00 


6 

1,31000 

-0.00 

0.00 

-0.00 

-0.00 

-12.24 

-450.00 


7 

.12500 

-0.00 

0.00 

-0.00 

-0.00 

-1.29 

115.16 

PRESTRESS 

ASSUMED FOR 

EIGENVALUE ANALYSIS 














PLATE 

NUMBER BREADTH 

NTV 

NLV 

NSV 

NTF 

NLF 

NSF 


1 

•12500 

-0.00 

3188.76 

-0.00 

-6.6o 

-6409.09 

334.84 


2 

2.58000 

-0.00 

3188.76 

-0.00 

-0.00 

-6409.09 

264.27 


3 

1,32800 

-0.00 

2816.46 

-0.00 

-0.00 

-5651.27 

185.73 


4 

1.25000 

-0.00 

13131.92 

-0.00 

-0.00 

-26268.11 

185.73 


5 

1,31000 

-0.00 

3188.76 

-0.00 

-0.00 

-6409.09 

450.00 


6 

1,31000 

-0.00 

3188.76 

-0.00 

-0.00 

-6409.09 

-450.00 


7 

.12500 

-0.00 

~28l6.46 

-0.00 

-0.00 

' -5651.27 

”115716” 


EIGENVALUE N O. LAMBD A _ FACTOR _ U IMIS ON f ACTOR 

1 m 1199999^ 6 EtOg . 23045 1309E»0l .230444102E?01 . .230451309E+0 1 

THE TOTAL AXI AL L OAD IS . 31330239E*05 . _ 

THE MODE IS _ .. _. _ - 

CLOCKWISE ROT. DEFLN. DOWN DEFLN. TO R T. L ONG . JDEF.L N , _ 

REAL PARTS ON ALTERN ATE LINE S* WITH IMAGINARY PARTS INSET BENEATH 

,9999997 IE- 01 -.2 59947 01 E ♦ 01 .8 3 8 1 1557E-03 .76889688E-0 1 

,00352 555 -.014 87365 ,0011 0046 .0U013180 

-.10000000E*00 -.259 ~ 94702E *01 -.8381 1443E-03 .76889692E-0 1 

• OQ 3525.53 ” ,0 1487369 • 001 1 0046 „-._p0013iei 

. 1000Q000E+0Q - ,259947 02E* 0 1 .838 1 1 443E-03 . 76889692E-0 1 

• 00352553 . 0 148 7 369 •00 1100 46 .000 131. 8 1 

- , 9999997 OE - 01 25994701 E + 0 1 -.8381 15 57 E - 03 • 7688 968 8E -01 

.00352555 .01487365 ~ .00il0046 -. 00013 180 

- ,59 1 031 2 OE - 1 1 -.268391 16E + 01 -• 17345568E-1 1 .75681 095E-01 

-,02001181 -.00000000 .00110677 0,00000000 

ALL REQUIRED EIGENVALUES HAVE BEEN FOUND FO R L AMBDA '« 32. 000000 

M LAMDA FACTOR "" 

1 .4000E + 01 .2388E + 01 _ 

? ,3200E*02 , 2304E* 0 1 


Table 5.- Input for Heated Titanium Panel 


24 36 7 200 500 2 

1 

BUCLASP-3 HEATED PANEL BUCKLING ANALYSIS 


1 1 

10000 

• 0 

2.0 

27.5 

-1. 

1. 

1. 

0 0 

0 2 

9 l 

0 2.001 

1.64E7 

1.64E7 

6.31E6 

• 3 

1.61E7 

1.61E7 

6.19E6 

• 3 

1.60E7 

1.60E7 

6.15E6 

• 3 

1.S6E7 

1.56E7 

6.00E6 

- .3 

1.52E7 

1.52E7 

5.85E6 

.3 

1.50E7 

1.50E7 

5.77E6 

.3 


1 .14 

2 *14 

3 .14 

4 .065 

5 .065 

6 .14 
-1 

1 

2 

3 

4 

5 

6 
-1 

1 1 U32 

2 2 *6575 

3 3 .48 

4 4 1,07 

5 5 1.07 

6 6 .9725 

9 2 .7025 

10 3 .4925 

-1 

4-74 270 

4 5 270 

512 1-9909000 

3 2 4 2 

3 3 4 3 

5 4 5 9 7 10 

3 5 6 1 

3 6-9909000 

3 7 8 7 

3 8 9 8 

3 9 10 6 


UUUUUUN* 


20 2 6 


ItO 


0.0 00000 

.0 

4.45E-6 

4.45E-6 

4.50E-6 

4.50E-6 

4.52E-6 

4.52E-6 

4.58E"*5 

4.58E-6 

4 . 65E-C 

4.65E-6 

4.69E-6 

4»69E-6 


80. 

80. 

100. 

100. 

175. 

175. 

250. 

250. 

280. 

280. 



Table 6.- Selected Output for Heated Titanium Panel 


THERMAL STRESS FOR AXIAL LOAD 

P = 0. 














NLF 

PLATE NUMBER BREADTH 

NTV 

NLV 

NSV 

NTF 

NLF 

NSF 

Cre £ +) 

1 

1*32000 

0.00 

0.00 

0.00 

0.00 

-1018.12 

0.00 

- 1020 

2 

• 65 750 

0.00 

0.00 

0.00 

0.00 

-186.06 

0.00 

-J&9 

3 

•46000 

0.00 

0.00 

0.00 

0.00 

19.19 

0.00 

IS 

4 

1.07000 

0.00 

0.00 

0.00 

0.00 

363.08 

0.00 

373 

5 

1.07000 

0.00 

0.00 

0.00 

0.00 

710.44 

0.00 

7i0 

6 

.97250 

0.00 

0.00 

0.00 

0.00 

1826.51 

0.00 

162D 

9 

*70250 

0.00 

0.00 

0.00 

0.00 

-188.06 

0.00 

- 169 

10 

.49250 

0.00 

0.00 

0.00 

0.00 

19.19 

0.00 

IQ 

PRESTRESS 

ASSUMED FOR EIGENVALUE ANALYSIS 






PLATE NUMBER BREADTH 

NTV 

NLV 

NSV 

NTF 

NLF 

NSF 


1 

1.32000 

0.00 

22960.00 

0.00 

0.00 

-45920.00 

0.00 


2 

.65750 

0.00 

22540.00 

0.00 

0.00 

-44268.56 

0.00 


3 

•46000 

o.co 

22400.00 

0.00 

0.00 

-43787.52 

0.00 


4 

1.07000 

o.oo 

10140.00 

0.00 

0.00 

-19467.28 

0.00 


5 

1.07000 

0.00 

9660.00 

0.00 

0.00 

-18611.45 

0.00 


6 

•97250 

0.00 

21000.00 

0.00 

0.00 

-39242.28 

0.00 


9 

.70250 

0.00 

22540.00 

0.00 

0.00 

-44268.56 

0.00 


10 

.49250 

0.00 

22400.00 

0.00 

0.00 

-43787.52 

0.00 


EIGENVALUE 

NO. lambda 

FACTOR 


LIMITS 

ON FACTOR 

TOTAL AXIAL 

1 

•27500 00 OOE +02 

.224138941 E +01 

•224133301E+01 

.224138941E+01 LOAD Cref.4 ) 

THE TOTAL AXIAL LOAD 

IS .44259060E+05 





.+5£15€*0S 


Reproduced from Reference 1 



Fig. 1 Basic load system for a plate (N^, and K g are forces per unit 



x 

Fig. 2 Defines the elastic properties (D^, D 12» D 13* D 22’ D 23* D 33» A ] 
A,_, A_„, A.,) of a plate, see also Fig. 1. 





45 


E 1 * 19*5xi0 6 

E 2 e 1 * 11x1 ° 6 ply thickness * 0.0052 

E 12 * 0.488x10° 

^ =0.19 



Figure 3.- Cross section of hat section panel used in example problems. 



PROBLEM 3C ( 0 -S - FLRNGES ) SHEfiR=45Q. 


Figure 4.- Machine plot of geometry for example problems. 



PROBLEM 3 C l 0 - S - FLANGES) SHLHR-45C). 
EIGENVALUE NUrl LAMBDA-4 -0000 


FAC TOR- 3 .7823 


7 

/ 


PROBLEM 3C i 0 - S - FLANGES ) SHEAR =4 50 . 
E I G E N VR L I J E N G - 1 LR M B UR - 4 . Q 0 0 0 



FRCT0R--3 .7823 IMAGINARY 



Figure 5.- Machine plot of buckling mode shape for example problem number 1. 

For problems with shear or anisotropy, the mode contains a real and an Imaginary 
part and both are shown. 


dUi-jt 


15 20 

LAMBDA 


PROBLEM 3C 'f 0 . 3 . PUHIC-.E6 ' SHLHR-4F.0 . 


EIGENVRLU 1 


Figure 6.- Machine plot of panel buckling strain as a function of buckle 
length for example problem number one. Strain ■ Factor xlO . 



